Preprint CHIBA-EP-105 

hep-th/9803133 
(revised) 
June 1998 



Existence of Confinement Phase in Quantum Electrodynamics 



Kei-Ichi Kondo 1 ^ 

1 Department of Physics, Faculty of Science, Chiba University, Chiba 263, Japan 
> E-mail: kondo@cuphd.nd.chiba-u.ac.jp 

Abstract 

We show that the four-dimensional U(l) gauge theory in the continuum formu- 
lation has a confining phase (exhibiting area law of the Wilson loop) in the strong 
coupling region above a critical coupling g c . This result is obtained by taking into ac- 
count topological non-trivial sectors in U(l) gauge theory. The derivation is based on 
the reformulation of gauge theory as a deformation of a topological quantum field the- 
ory and subsequent dimensional reduction of the D-dimensional topological quantum 
field theory to the (D-2)-dimensional nonlinear sigma model. The topological quan- 
tum field theory part of the four- dimensional U(l) gauge theory is exactly equivalent 
to the two-dimensional 0(2) nonlinear sigma model. The confining (resp. Coulomb) 
phase of U(l) gauge theory corresponds to the high (resp. low) temperature phase of 
0(2) nonlinear sigma model and the critical point g c is determined by the Berezinskii- 
Kosterlitz-Thouless phase transition temperature. The quark (charge) confinement 
in the strong coupling phase is caused by the vortex condensation. Thus the con- 
tinuum gauge theory has the direct correspondence to the compact formulation of 
lattice gauge theory. 
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1 Introduction 



In this paper we study the phase structure of the continuum Abelian U(l) gauge 
theory by including the effect due to compactness of U(l) group. The reason of taking 
into compactness is as follows. From the viewpoint of unified field theory, the Abelian 
group should be embedded as a subgroup in larger non- Abelian gauge group. In view 
of this, the Abelian group should be compact. Another importance of compactness 
of Abelian gauge group stems from a possibility of explaining quantization of charge 
|]T|. In non-compact QED there are no reasons for charge quantization. 

In this paper we show that the four-dimensional U(l) gauge theory has a confine- 
ment phase in the strong coupling region g > g c due to the compactness (periodicity) 
leading to the non-trivial topological configuration. If we neglect the periodicity, 
we have a free U(l) gauge theory which has only one phase, the Coulomb phase, 
as expected. This work confirms the claim made by Polyakov |2|, |3|. However, the 
claim that the abelian gauge theory has a confinement phase sounds strange from the 
conventional wisdom based on the continuum abelian gauge theory. We clarify the 
meaning of this statement in what follows. 

More than twenty years ago, it has been pointed out by many authors that the 
four- dimensional SU(2) non- Abelian gauge theory bears many similarities with a two- 
dimensional 0(3) nonlinear sigma model (NLSM). Both theories possess asymptotic 
freedom, multi-instanton (and anti-instanton) solution, dynamical mass generation 
and scale invariance (i.e., no intrinsic scale parameter), see H for references. 

These similarities can be seen also in the lattice regularized versions of these 
models, between spin models and lattice gauge theories ||. Naively the scaling limit 
of the classical 0(3) Heisenberg model is the 0(3) NLSM, whereas that of the SU(2) 
lattice gauge theory is the SU(2) gauge theory. One can take the scaling limit of a 
lattice theory at a second order phase transition point. Hence the scaling limit is 
taken by approaching the critical point T — > T c (or g — > g c ) as the lattice spacing a 
goes to zero, a — > 0, in such a way that the physical quantities remain finite. 

In two-dimensional classical 0(3) Heisenberg model, the two-point correlation 
function decays exponentially at any finite temperature. This corresponds to a claim 
in four-dimensional lattice SU(2) gauge theory that confinement phase survives as 
long as the coupling constant g is positive, even if g <C 1. Both models have a phase 
transition at T = (g — 0), i.e., T c = (g c = 0) which is believed to be the second 
order. 

In a previous paper M, it has been shown that these similarities between two 
models are not merely an accident; actually we have proved the exact equivalence 
between (D-2)-dimensional 0(3) NLSM and the D-dimensional topological quantum 
field theory (TQFT) obtained by removing the perturbative deformation (topological 
trivial sector) from the D-dimensional SU(2) non- Abelian gauge theory (D > 3). 
This proof is based on an idea of the dimensional reduction of Parisi and Sourlas ||. 
The case of D = 4 is the most interesting case of physical reality. 

What can we say in the Abelian case? For this, recall a fact that the two- 
dimensional 0(2) NLSM or XY model undergoes a phase transition without the 
appearance of a spontaneous magnetization. This absence of order parameter in two 
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dimensions is consistent with the Coleman-Mermin- Wagner (CMW) theorem 0. The 
low temperature phase (T < T c ) contains massless spin waves. On the other hand, the 
high temperature phase (T > T c ) is completely disordered. For this phase transition, 
the periodicity of the angular variable if is quite essential. The model has topological 
singularities, called the vortices. These vortices condense in high temperature and 
disorder the correlation function ||. This phase transition is called the Berezinskii- 
Kosterlitz-Thouless (BKT) transition || . The vortex part is equivalent to the neutral 
Coulomb gas and sine-Gordon theory, see 0, [10, 11]]. Although the existence of BKT 



transition is rather subtle, the existence of BKT transition was rigorously proved by 
Frohlich and Spencer ||12|| . 

In lattice formulation, it is well known that all of these properties in two-dimensional 
abelian spin models have correspondences in the Abelian gauge theory in four dimen- 
sions. The vortices in two dimensions are closely related to magnetic monopoles 
in four dimensions, see [13, T3, |T5| . The condensation of closed loops of magnetic 



monopole currents leads to quark (charge) confinement in the strongly coupled phase 
of Abelian gauge theory, since electric flux cannot easily penetrate such a medium 
(which we call the dual Meissner effect). It will be worth remarking that the dual su- 
perconductor vacuum of quantum chromodynamics (QCD) has been derived recently 
without any ad hoc assumption [16] from QCD in the continuum formulation. 

In lattice gauge theory, charge confinement in the sense of area decay of the Wilson 
loop is derived in the strong coupling region by using the strong coupling expansion 
H [17|. Quite remarkably, the quark (charge) confinement in lattice gauge theory 
occurs irrespective of the details of the gauge group, as long as it is compact (discrete 
[|18| , or continuous), even for abelian gauge group. However, one expects that the 
U(l) lattice gauge theory in four dimensions (£7(1) 4 ) has a Coulomb phase in the 
weak coupling region, which was proved rigorously by Guth |2(| and Frohlich and 



Spencer plR . Therefore the U(l) lattice gauge theory undergoes a phase transition 
at a finite non-zero coupling g c . In continuum gauge theory, such non-trivial phase 
structure was suggested to occur due to topological non-trivial configurations by 
Polyakov P, |TJ. Actually he has shown the confinement phase in three-dimensional 
U(l) gauge theory for arbitrary gauge coupling, in agreement with the lattice analysis. 
In four dimensions, he claimed that the weak coupling U(l) gauge theory does not 
confine. Accordingly, it is expected that the U(l) gauge theory in four dimensions 
have two phases, confinement and deconfinement (Coulomb) phases, whereas only 
one phase i.e., confinement phase exists in three dimensions. 

In this paper, we show that the continuum four- dimensional U(l) gauge theory has 
two phases, (strong coupling) confinement phase and (weak coupling) Coulomb phase, 
which have direct correspondence with the high-temperature and low-temperature 
phases in 0(2) NLSM respectively. The phase transition point corresponds to the 
BKT transition in XY model. Therefore, the phase transition point g c is determined 
by the BKT transition temperature T c . This is one of the main results of this paper. 
This result is obtained as a specific case of the previous paper [[|. 

Therefore, in the strong coupling phase (g > g c ), the continuum U(l) gauge theory 
confines quarks and the gauge field becomes massive, in agreement with the result 
of lattice gauge theory. In the weak coupling phase, on the other hand, quarks are 
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liberated and gauge field remains massless. In the weak coupling phase (g < g c ), the 
beta function of the renormalization group is identically zero and < g < g c is 
the line of fixed points, if we neglect the perturbative deformations. For a review of 
lattice gauge theory, see e.g. |23|, |24|, |25| . 



The plan of this paper is as follows. In section 2, we give the reformulation of 
Abelian gauge theory as a deformation of a topological quantum field theory. In sec- 
tion 3, using the reformulation of section 2, we evaluate the Wilson loop expectation 
value in four-dimensional Abelian gauge theory. In the final section we discuss the 
renormalization group properties and extension of our scheme to other dimensional 
cases. More details about the interplay between the Abelian and non- Abelian cases 
are given in a forthcoming paper. 



2 Abelian gauge theory as a deformation of TQFT 
and dimensional reduction 

Now we reformulate quantum electrodynamics (QED) as a deformation of a topologi- 
cal quantum field theory. It is obtained as a special case of non- Abelian gauge theory 
given in the previous paper [|]. 



2.1 Decomposition into perturbative and topological non- 
trivial sectors 

QED on the .D-dimensional space-time is defined by the action, 

Sqed = Jd D x(C QED [a,^]+C GF ), (2.1) 



£qsd[cv,V] := —^UuUu +^{iYDM -m)ip, (2.2) 



where 



ftj,u{x) := df.a^x) - dvd^x), (2.3) 

DM := d^-iga^. (2.4) 

The gauge transformation of the U(l) gauge field a M (x) and the fermion field ip is 
defined by 

a^x) -> aV( x ):=a^x) + -U{x)d^U\x), U (x) G U (1), (2.5) 

ip(x) ip u (x) := U(x)ip(x). (2.6) 
The gauge fixing term C GF is given by 

£gf ■= -i8 B G gf [a^ C, C, 0], (2.7) 



3 



using the nilpotent Becchi-Rouet-Stora-Tyupin (BRST) transformation 5b, 

daa^ix) = <V7(:z), 

5 B C(x) = 0, 

5bC(x) = i(j)(x), 

5 B 4>{x) = 0, 

5 B ip(x) = igC(x)ip(x), 5 B ${x) = -igC(x)ip(x), (2.8) 

where <ft is the Lagrange multiplier field. 

The partition function of QED with the source term 

Sj[a^ C, C, 0, V, j>\ ■= J d D x(J fl a^ + J C C + J- C C + + # + rrf) (2.9) 

is given by 

Zqed[J] ■= J [da?] [dC] [dC] [d<j)\ [#] [#] exp [iS% ED + iSj) . (2.10) 

To reformulate QED as a deformation of a topological quantum field theory ac- 
cording to |Q, we first regard the U(l) gauge field a M and the fermion field ip as the 
gauge transformation of the U(l) gauge fields and 

aJx) : = vjx) +uJx), u)Jx) := -U{x)d l JJ^{x). (2.11) 

9 

ip(x) := U(x)^(x), (2.12) 

Here [] and \& are identified with the field variables in the perturbative (topological 
trivial) sector (Q — 0), whereas belongs to the topological nontrivial sector (Q ^ 
0). 

Furthermore we introduce new ghost field 7, anti-ghost field 7 and the Lagrange 
multiplier field j3 in the perturbative sector. They are subject to a new BRST trans- 
formation 5b, 

8 B l{x) = 0, 
Sbj(x) = i(3{x), 
5 B (3(x) = 0, 

5 B ^(x) = igj(x)iff(x), 5 B ^(x) = -igj(x)^(x), (2.13) 
Then the partition function of QED is rewritten as 

Zqed [J] = J [dU] [dC] [dC] [d(j>] J [tfyj [dj] [dj] [d(3] [dV] [dV] 
x exp |i J d D x( - iSsGgflu^ + u M , C, (7, 0]) 

+z /" d D x(£ QSD [w, - i5 B G gf ^,7,1,0]) 
+iSj[u^ + v„ C, C, cf>, UV, $ £/t]}, (2.14) 

1 The decomposition of a„, a„ = + a;„, corresponds to the superposition of two independent 
configuration, (p = ipsw + Vv (spin waves and vortex parts) in XY model. 
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where G g f is a gauge fixing functional for the perturbative (topological trivial) sector. 

2.2 Gauge fixing 

The Lorentz gauge is given by 

F[a] := = 0. (2.15) 
The most familiar choice of G g f 

G gf = C(dX + f 0) (2-16) 
yields the familiar form of the gauge fixing term, 

— — Ol 

Cgf ■= -i6 B G gf [a^C,C,(l)} = (l)dX + iCd> i dt l C + -(l) 2 . (2.17) 
In this paper we propose to use a choice, 

<%W = -h(\<i + iCCy (2.18) 
where S B is the anti-BRST transformation H, 



5 B a^(x) = M C(ar), 

S B C(x) = i<f>(x), 

5 B C(x) = 0, 

5 B (j)(x) = 0, 

5 B ip(x) = C(x)ip(x), 

(f){x) + (j)(x) = 0, (2.19) 

where is defined in the last equation. 

Apart from a total derivative term, this choice yields 

C GF = iSjsh QaJ + iCc\ = -i5 B [C{d^ - 0)]. (2.20) 



Therefore, the choice ( 2.18Q corresponds in ( |2.17| ) to the choice of the gauge-fixing 
parameter, 



a = -2, (2.21) 

I also in the 

the decomposition, 



which has appeared also in the non-Abelian case Q. The above choice for G^j 1 ^ yields 



Car = -iSBG^l^ + v^CCA] (2.22) 
C-tqft + iv^5 B 5 B u^ (2.23) 



iS B 5 B ( -(up + t^) 2 + iCC 



where we have defined 

1 



Ctqft ■= i5 B 5 B [^l + iCCy (2.24) 

Here we have used that the action of 5 B is trivial in the perturbative sector, 

S BVli = = 5 B Vp, (2.25) 

while 

5 b uj^ = d^C, 5 b uj^ = d^C. (2.26) 

2.3 Deformation of topological quantum field theory 

Finally, the partition function of QED is cast into the form, 

Z Q ed[J] ■= J[dU][dC\[dC\[d^ V {iS TQFT [u^C,CA\ 

+i J d D x[ru fl + j c c + j- c c + j^>] + iw[U; r,v,r]}}, (2.27) 

where W[U; J^,fj : r]\ is the generating functional of QED in the perturbative sector 
(pQED) given by 

e iw[u;j^,v] := j [^][d 7 ][d7][^][^][^]exp{^ pQED [^,*, 7 ,7,/3] 

+1 J d D x[v^ + fjUV + riVU*]}, (2.28) 

SpQEDfy,*, -y,Tf,0\ ■= j d D x[C QED [vM -iSsGgfiv^-y, 7, f3)], (2.29) 
J„ := J^ + iSbSbU^. (2.30) 

The correlation functions of the original (fundamental) field a^,^,^ are obtained by 
differentiating Zqed[J] with respect to the corresponding source ^,77,77. 

All the field configurations are classified according to the integer- valued topological 
charge Q which is specified later. The above reformulation of gauge theory is the 
decomposition of the original theory into the topological trivial sector with Q = 
and topological non-trivial sector with Q ^ 0. This corresponds to the decomposition 
of XY model into spin wave part (Q = 0) and the vortex part (Q 7^ 0) where Q is 
given by the winding number of the vortex solution. However, the XY model is not 
a gauge theory and does not have any local gauge invariance. 

The integration over the fields (U, C, C, <fi) in TQFT should be treated non- 
perturbatively by taking into account the topological non-trivial configurations. The 
deformation W[U; J M , fj, rj] from the TQFT should be calculated according to the or- 
dinary perturbation theory in the coupling constant g. The perturbative expansion 
around the TQFT means the integration over the new fields (v^, 7, 7, (3) based on the 
perturbative expansion in powers of the coupling constant g. 
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2.4 Dimensional reduction to 0(2) NLSM 

Following the argument given in M based on the Parisi-Sourlas dimensional reduction, 
it turns out that the D-dimensional TQFT (as the topological non-trivial sector of 
D-dimensional U(l) abelian gauge theory) with an action 

Stqft[oo^ C, C,(f>] = J d D x i5 B 6 B (^u^u^x) + iC(x)C(x)^j (2.31) 
is equivalent to the (D-2)-dimensional 0(2) NLSM with an action, 

So(2)nlsm[U] := 2tt J d D - 2 z ^(z)uj^z), u^z) := -U{z)d^U\z), (2.32) 

= / d D - 2 z ^d,U(z)d,U\z). (2.33) 
J 9 



Dimensional reduction is due to a fact that the action ( p. 31 ) has a hidden super- 



symmetry and can be rewritten in the 0Sp(D/2) symmetric form in the superspace 
formulation, see M. 



3 Quark confinement in Abelian gauge theory 

Now we calculate the Wilson loop expectation in the U(l) gauge theory based on the 
reformulation given in the previous section. In what follows we move to the Euclidean 
formulation. 

3.1 Dimensional reduction of Wilson loop 

We define the Wilson loop operator for the closed loop C by 

Wc[a] = exp (iq j> a fl (x)dx fl \ , (3.1) 
where q is a test charge. In the abelian gauge theory, the Wilson loop factorizes, 

Wc[a] = exp (iq <j> uj^(x)dx^j exp (iq j> v )1 (x)dx^j =: W c [a;]WcM- (3.2) 
For Q = sector, we choose the gauge-fixing function, 

Arf = 7(0X+|9), (3.3) 

with a gauge-fixing parameter £. For the U(l) gauge theory with an action (omitting 
matter fields), 

Spu{i)[v,j,j,0\ ■= -^(d^v - d^v u f - iS B G gf , (3.4) 
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the perturbative part is given by 

e lW[uJ > JA01 = J [tfyj [dj] [di\ W] exp { - V(D [v, J,JiP}+iJ d D xv,J,). (3.5) 
Integrating out the fields 7, 7, /? yields 

,iWMo,o] _ /" exp | _ s[v] _ f d D xtv fl (x)[J ll (x) + i5 B 5 B ^(x)]\, (3.6) 



e 



where 



S[v] := / d D x 



(3.7) 



For the calculation of the Wilson loop expectation, J M is taken to be the current 
along the closed loop C such that 



d D xv fJ ,(x)J^(x) = q j> v fl (x)dx^ 1 . (3i 
It is easy to see that 

giWKJ.0,0] = (WcW^ Sb ' Sb ^) pU{1) J [dv^]e- s[v]+{v ^ B ~ &B ^\ (3.9) 

where 

{vh,6 b 6bu)h) ■= J d D xv fl (x)S B S B ^ fl {x). (3.10) 
The Wilson loop expectation is rewritten as 



(W c {Lu}e lW ^> J > ^) 



where 

e tW[ w ,o,o,o] = ( e ^' SB ~ SB ^) pm f[dv fl ]e- s ^ +( -^' 5B ' SB ^. (3.12) 
Expanding the exponential 

e (v^,S B SB^n) = e J d D xS B 5 B (vnUn)( x ) _ e J d D x{Q B ,d B (v ll ui li )(x)} /g 

into power series and using a fact that the vacuum of TQFT obeys 

Qb\0)tqft = 0, (3.14) 



we find that this term does not contribute to the expectation value ( p. 11]) . Therefore 



the Wilson loop expectation is completely separated into the topological part (Q ^ 0) 
and the perturbative part (Q — 0), 

(W c [a}) u{1) = (W c [uj}(Wc[v}) pU{1) )tqft = {W c [uj})tqft(Wc[v}) pU{1) . (3.15) 



S 



This corresponds to the result eq.(10) of Polyakov 0. This property does not hold 
in the non-Abelian case, which makes the systematic calculation rather difficult. 

In order to use the dimensional reduction for calculating the topological part, 
we choose the loop C so that C is contained in the D — 2 dimensional space. For 
D — 4, the loop C must be planar. Then the dimensional reduction of the topological 
part leads to the equivalence of the Wilson loop expectation value between the D- 
dimensional U(l) TQFT and (D-2)-dimensional 0(2) NLSM, 

(W c [uj{z)]) T q FTd = (Wc[co(z)])o(2)nlsm d . 2 , z G R D ~ 2 . (3.16) 

if the Wilson loop has its support on the (D-2)-dimensional space on which the NLSM 
is defined. Hence, the calculation of the Wilson loop in the four-dimensional U(l) 
gauge theory is reduced to those in the two-dimensional 0(2) NLSM and the four- 
dimensional perturbative U(l) gauge theory, 

(W c [a}}u(i) D = (Wc[u])o(2)nlsm d - 2 (Wc[v]) p u(i) d - (3.17) 
For large rectangular loop with sides R, T, the static potential is obtained by 

V(R)=\\m ^\n(W c [a}). (3.18) 

1 — >oo _/ 



In four dimensions (D — 4), it is well known p3f that for large Wilson loop (Wc[v]) p u(i) 4 
gives the Coulomb potential, 

q 2 1 

V(R) = - — — + constant. (3.19) 

47T R 



For a derivation, see e.g. Appendix of pq| . 

In the following, we show that (Wc[uj})o(2)nlsm 2 exhibits area law for strong cou- 
pling g > g c with a finite and non-zero value of a critical point g c . This confinement- 
deconfinement transition corresponds exactly to the BKT transition. 

3.2 0(2) NLSM and Wilson loop 

Defining the angle variable (p(z) for U(z) G U(l), 

U(z) = e^ z \ (3.20) 

we obtain 

u,(z) = -U{z)d,U\z) = -d^{z). (3.21) 
9 9 

Then the action of 0(2) NLSM reads 

]D-2 /_\, . /_\ P f jD-2, o / \o / \ n 27T 



So(2)[U] = / d u - 2 zixuj,(z)oJ^z) = ^ / d u - l zd^{z)d^{z), (3 := — . (3.22) 
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The partition function is defined by 

Zo(2)[U] = I [dU] exp(-S 0(2) [C/]), W] = II (3-23) 

using the Haar measure dU on U(l). For the notation of the two-dimensional vector, 

S(x) — (cos(p(x), sixi(p(x)), (3.24) 

the action reads 

So(2)[U] = £j d D - 2 z <9 M S(x) ■ d,S(x). (3.25) 

Hence the 0(2) NLSM is regarded as a continuum version of classical planar spin 
model. In the following, we identify (3 with the inverse temperature 1/T. Hence the 
high (resp. low) temperature of the spin model corresponds to strong (resp. weak) 
coupling of the gauge theory. 

It should be remarked that 0(2) NLSM with the action ( |3.22| ) is not a free scalar 
field theory, since this theory is periodic in the angle variable <p (modulo 2tt). Of 
course, if we neglect this periodicity and treat the variable ip as a non-compact variable 
<p(x) G (— oo, +oo), we have a trivial theory, i.e., free massless scalar field theory. In 
this case, the Contour integral is zero, 

u^{z)dz^ = -f dpcpWdz 1 * = 0. (3.26) 

Hence, the Wilson loop W c [u] is trivial, W c [w] = 1 and hence the total static quark 
potential comes from the the Wilson loop expectation W c [v] of perturbative U(l) 
gauge theory and is equal to the Coulomb potential. Thus we obtain a trivial result 
that the four-dimensional non-compact Abelian gauge theory fails to confine quarks 
(charges). 

However, the periodicity (or compactness) leads to topological non-trivial solu- 
tions which are seeds for confinement, as shown in what follows. 

In the following, we restrict our consideration to d := D — 2 = 2 case. The 
extremum of the classical action (|3.22| ) is obtained as a solution of the classical field 



equation, 

VV = (mod 2?r). (3.27) 

The harmonic function (p is constant or has singularities. We require (p be constant 
at infinity and assume only isolated singularities around which (p varies by ±2ir as 
one turns anticlockwise. Then the solution Q is written as 

t z _ z .\ 

(p(z) = Qi arctan l — = Qilm\n(z — zA, z:=Xi+ix 2 - (3.28) 

i (z-Zi)i i 



2 The two-dimensional Laplace equation is equivalent to the Cauchy-Riemann equation. Hence 
the solution is given by the holomorphic function. If one avoids branch cuts, it is a meromorphic 
function. 
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This denotes a sum of vortex excitations located at points x% G R 2 and of vorticity 
Qi (integers). The solution has an alternative form, 



¥>(*) = In *t! =E 



In ( z Z ^ + ) _ ln ( z ~ Z A ) 



(3.29) 



l z - z< I |*~ *i 

This means that vortices of intensity ±1 are centered at the points zf and that 
intensities of higher magnitude are obtained when several zf (or z~) coincide. Note 
that the angle <p is a multi- valued function, but e %Lp is well-defined everywhere, except 
at the singular points. 

The contribution of the solution (|3.28|) to uj^ is 

Un{z) = -d^{z) = - ^2 QiE^ Z _ Zl " = - ^Qit^ v d v \\i\z - Zi\. (3.30) 
9 9 i \ z z i) 9 i 

Therefore, the integral of one- form u := u^dx^ along the closed loop C is 

u^da? = I u = - r dQi = £ —Qi, (3.31) 
c Jc g t Jo 7 9 

where the sum runs over all the vortices inside the closed loop C and B« is an angle 
around z = z iy 

Qi(z) : = arctan ~ Zi \ 2 . (3.32) 

[Z — Zi)\ 

Note that u is a closed form, du = 0, but it is not an exact form, that is, there does 
not exist a function (zero-form) such that u = df with / being defined everywhere in 
R 2 — {0}. Domain of f(z) = Q(z) is restricted to R 2 — R + , in other words, for one 
unit vortex at the origin, 

x 

= e ^ v ^I ~ 2vr6»(xi)(5(x2)^2- (3.33) 

This is analogous to the case of the magnetic monopole in three dimensions where 
the magnetic field is given by 

H » = \r j k ~ 2 ^sA^(x 2 )e(x 3 ). (3.34) 

The singular line of ( |3.33|) in two dimensions does not contribute to the action, so 
does the Dirac string (on the positive Z axis) in three dimensions. 

In order to calculate the classical action for the singular configuration ( |3.28| ), we 
consider a disk of radius Rq, Di := {\z — zi\ < R ; z 6 R 2 } (centered on each singular 
point Zi) which is small with respect to the distances between vortices. Let 71 be the 
remaining domain of integration outside the vortices. The classical action consists of 
two parts, the self-energy (action) part of vortices, 

S« = §W d 2 z(V V (z)) 2 , (3.35) 

Z ■ J\z-Zi\<R 



11 



and the remaining part, 

P 



S (2) := § / (PziV^z)) 
2 Jn 



= -2n[3 Q l Q J In \ Zi - Zj \ +J2 Qi*P ln l / R o- (3-36) 

i¥=3 » 

Summing over all vortex sectors leads to the partition function of the form, 



oo £n - n 
n=0 J j=l 



(27r) 2 /3^g^A(^,^) 



ho 



C:=e~ s(1) , (3.37) 



where C comes from the self-energy (action) part of vortices and A expresses the 
two-dimensional inverse Laplacian given by 

A(x,0) = — ln-^-. (3.38) 
2n \x\ 

Therefore the partition function just agrees with the two-dimensional neutral Coulomb 
gas (i.e., a gas of classical charged particles with a Coulomb interaction and globally 
neutral, J2i Q% — 0)- 

The transition temperature is estimated as follows. The contribution to the free 
energy from one vortex pair at distance r 12 in a box of linear dimension L is 

F ~ ln / d 2 z 1 d 2 z 2 exp(-2irl3\n\z 1 - z 2 \/R ) 

J\z 1 -z 2 \>Ro 

~ ln[L 4 exp(-27T/31nL/ J R )] 

~ (4-27r/3)lnL. (3.39) 

The vortices always arise in pairs of opposite Coulomb charges to yield finite energy 
configurations and each pair forms an elementary dipole (Coulomb dipole gas). If 
2-7T/3 > 4 (the low temperature or weak coupling phase), the contribution from the 
vortex pair is negligible in the limit L — > oo. In this phase, charges are bound and one 
has a dielectric medium. As (3 — > oo, few vortices are present and their correlation 
decreases rapidly with the relative distance. This describes a dielectric medium of 
neutral bound states. In this regime, the correlation (U (z)U(z')) decays polynomially, 

(U(z)U(z')} = \z- z'\-**. (3.40) 

On the other hand, if 2-7T/3 < 4 (the high temperature or strong coupling phase), an 
instability occurs and the creation of well-separated vortices is favored and disorder 
increases. In high temperature phase, one has a plasma of free charges. The vortex 
expectation decays exponentially yielding exponential decay of (U(z)U(z')), 

(U{z)U{z')) = \z- z '\-ihe- m ^ z - z 'l (3.41) 
Hence a naive estimate of the critical temperature is obtained 

& = -, 9 2 e =* 2 . (3-42) 

71 
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This is the phase transition without the appearance of a spontaneous magnetization. 
The phase transition can be interpreted as a dipole condensation. The critical point 
separates the dissociated dipole phase from the condensed phase. 

In order to calculate the Wilson loop, we use the equivalence of the Coulomb gas 
to the sine-Gordon model (see Appendix A for a proof). The sine-Gordon model is 
defined by the action and the partition function, 



S s g(4>) '■= oc d d x 



^(^0(x)) 2 - hcos(p{x) 



(3.43) 



Z sG {h) := J[d(j)\e-x V [-S sG {(j>)}. (3.44) 

This is equivalent to the partition function of a globally neutral gas of particles of 
charges Qi = ±1 through a Coulomb potential in d dimensions, 

oo 2n n , i 

z c ,±i(h) = Eri nMV%ex P {-JE[%-^l)+%.-!/ 3 l)] 

n=0 \ n -) i=l 1 i<j 

-£V(I*<-Vil)]} (3.45) 
at temperature T with the fugacity z, 

a = T= ik = ^' " h = 2 * = *- (3.46) 

Note that a/3 = l/(4ir 2 ). It is known that the transition point of sine-Gordon model 
is 

ac =±, (3.47) 



which is in agreement with (|3.42 ). The relation of tp and <f> is given by 

/■OO 

<f)(x,t) = / dy<p(y,t), (3.48) 

or 

= e^cj). (3.49) 

This implies that the fields (p and <ft are dual variables. Therefore, 0(2) NLSM 
is equivalent to the Coulomb gas and moreover it is equivalent to the sine-Gordon 
model when charge of the Coulomb gas (or vorticity of 0(2) NLSM) is restricted to 
Qi = ±1. Taking into account Qi > 1 will lead to the cos(Q0) term. The above 
consideration can be transferred into the lattice formulation, see 112 311. 



3.3 Wilson loop and area decay 

The Wilson loop expectation is calculated using the equivalent sine-Gordon model. 
The generating functional for the charge density p 

p(x) := ^QiS^ix-Xt), (3.50) 
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is obtained as 



141 = ^i^^) sG =^^) sG , 



Z sG [rj] = J[d(f)}exp 



a Id x 



(3.51) 

\(d^(x)) 2 - h cos[0(x) + J7(x)]] } . (3.52) 



In two dimensions, we can introduce the dual vector field, 

Hfj, = e^Uy. 

Then the dual field is connected with the charge density p as follows, 

9 

Note that the rotation of uo^ or the divergence of 



(3.53) 



£u^z)dz^ = j^d^z^z = J d II H II (z)d 2 z = — J p(z)d 2 z. (3.54) 



(3.55) 



measures the density of topological charge. If we identify the right-hand-side with 
the magnetic charge, this implies Dirac quantization condition, 



27T 

g m = —Q (Q ■ integer). 

The Wilson loop is calculated from the i] given by 

[z - x) v 
[z — x) 2 ' 



q f i 
r](x) = - f dz^e^u- 
9 Jc 



since 



q <j)^u fl dz fl = J d d xp(x)r](x) = Y^QiV( x d- 



(3.56) 



(3.57) 



(3.58) 



Note that r)(x) = if the argument x of r)(x) is outside the loop, while rj(x) = 2nq/g 
if x is inside the loop. 

In high temperature phase, photon is massive, whereas photon is massless in the 
low temperature phase. This is because in the high temperature phase, random distri- 
bution of free vortices with long range interaction spoil the correlation. Therefore, in 
the high temperature phase, the Wilson loop expectation is estimated by the steepest 
descent as 

(WcM) sG = exp|-a J d 2 x[^(^[0 d (x)-77(x)]) 2 -/icos0 d (x)]J, (3.59) 

where (f> c i is determined by the Debye equation, 

V 2 [Mx) - v(x)] = hsm<p cl (x). (3.60) 
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Corrections to this field due to fluctuation are exponentially small in high temperature 
phase. The loop C is placed in two-dimensional plane. We can perform the calculation 
in the same way as done by Polyakov |3|]. 

Instead of repeating similar calculation to we use the Villain form [[2~7jj , 



e Jcos0 _^ e J p -|(^-27rm) 2 



]T e~^- 2 ^ , (3.61) 



tn6Z 



to estimate the Wilson loop expectation. Then the partition function is replaced with 
(apart from field-independent constants) 



Z sG [ V ] = /[d0]e- a / d M§(^)) 2 ] rj e 

= f[dcl>}e- a f ddx l^<f>(*)) 2 ) Yl J2 



ah cos[<f)(x)-\-ri(x)] 

&R d 

[<t>(x)+r)(x)-2Trm(x)] 2 

x&R d m(x)£Z 

J2 [[d(j)}e~ a f d d x\[-^{x)d 2 <j>{x)} e -^ I ^xi^+^x^Mx)? 

{m(x)£Z;x£R d } 

S d d x[ri(x)-2nm(x)} 2 

{m(x)£Z;x£R d } 

X J[d(j)}e~ I d d x{f4>(x)(-d 2 +h)cf>(x)+o l h<l>(x)l V (x)-2nm(x)]} 

exp — — / d d x{[rj{x) — 2nm(x)] 2 

{m(x)eZ;xeR d } ' 

-h[r)(x) - 2vrm(x)](-<9 2 + h)- l [r](x) - 2nm(x)]}}. (3.62) 



When f] = 0, the denominator is obtained, 



Z sg[0] = Yl ex P 

{m(x)£Z;x£R d } 



-(2vr) 2 ^ I d d x{m{xf - hm(x)(-d 2 + hy l m(x)} [3.63) 



Note that the field 77(0;) has its support on S(dS = C) and has the value 2n(q/g). 
If q is an integral multiple of g (the elementary charge), we have 77 G 2nZ. This 
is absorbed by the shift of m. Therefore, in this case, charge confinement does not 
occur. This is interpreted as the charge screening. 

A naive estimate of the ratio Z s g[i]]/Z s g[0] is given when rj ^ 2irZ in Appendix 
B. Finally we obtain the area decay of the Wilson loop expectation, 

(WcM) sG = e~^ c \ (3.64) 

This implies the linear static potential 

V(R) = aR (3.66) 
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9c 




Figure 1: Renormalization group beta function of U(l) gauge theory 



between two fixed electric charges and electric string with uniform energy density a 
which is called the string tension. Therefore condensation of topological non-trivial 
configuration leads to quark confinement. From ( |3.17| ), ( |3.19| ) and the total 

static potential is given by 

q 2 1 

V(R) =aR- — — + constant. (3.67) 

Even the continuum Abelian U(l) gauge theory has a confinement phase, exhibit- 
ing rich phase structure as in lattice compact U(l) gauge theory. 



4 Discussion 

In the four-dimensional pure U(l) gauge theory, we have proved the existence of a 
strong coupling phase where the fractional electric charge is confined by the linear 
static potential due to vortex condensation. In the following we discuss a few points 
of perspective. 

4.1 Renormalization group and non-Gaussian fixed point 

The O(N) NLSM has the renormalization group beta function, 

P{g) : = /jM = .^ + 0(g % (4.68) 

In the low temperature (0 < T < T c ), therefore, the 0(2) NLSM has vanishing 
beta function, (3 = (T < 1), and < T < T c is the line of fixed points. This is 
consistent with the fact that at low temperature, the inverse correlation length or mass 
m = -T 1 of 0(2) NLSM or XY model vanishes, i.e., m{T) = for all T < T c . The 
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theory is conformal invariant. For high temperature (T > T c ), the renormalization 
group study of XY model shows that the mass behaves as 

mCO~exp(- ^r^ ), TIT C , (4.69) 

with a constant C. 

These results established in two-dimensional model would be translated into the 
four-dimensional Abelian gauge theory, provided that two theories have the same 
renormalization-group beta function. We assume that in two theories the mass m(g) 
is generated by the dimensional transmutation in such a way that the beta function 
(3(g) is related to the mass m(g) through the well-known relation, 

m(g) = fif(g) = fi exp (- f , (4.70) 

where (3(g) is defined by 

dg(fj) f(g) 

^-=^--wr (471) 

In the weak coupling phase (g < g c ), the gauge field is massless, m(g) = and the 



beta function of the renormalization group [22] is identically zero 



(3(g) = (0<g<g c ). (4.72) 

Therefore, < g < g c is the line of fixed points. In strong coupling phase (g > g c ), 
the beta function behaves as (see Figure 1) 

M = ~^-(9 2 ~ & 2 < 0, {g> g c ,g = g c ), (4.73) 



which is compatible with the BKT mass ( |4.69j ) 



m(g) = Aexp\ I , g j g c , (4.74) 



9 2 ~9c 



with constants A, B > 0. Note that (3(g) is independent of A. 

It is worth remarking that the recent lattice computer simulation [^, ^| indicates 
the existence of an non-Gaussian fixed point in four- dimensional pure compact U(l) 
gauge theory. This should be compared with the old results BO, |31l| . In the simulation 



[p8j p29f the continuous phase transition was found and analyzed according to the 
power law scaling, although our investigation suggests the scaling behavior of the 
essential singularity type and the data do not exclude the essential singularity. It will 
be rather difficult to specify the essential singularity in the computer simulations, 
since the lattice size available is not yet so large to confirm this issue. Anyway, it will 
be interesting to find any relationship to fill the gap between two approaches. 
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If the mass scale is generated by dimensional transmutation, the above results are 
quite analogous to the situation found for the dynamical fermion mass and the beta 
function in quenched massless QED The relationship of these results with the 

quenched QED is more suggestive by the method of bosonization or fermionization. 
The two-dimensional sine-Gordon model is equivalent to the massive Thirring model 



34, RHI with an action, 



Sty,® = J d 2 a#(*7A + mty - j(^) 2 ]. (4.75) 

The correspondence between two theories is given by 

G 

l + _ = Ana, (4.76) 

7T 

$1^ = —^e»vd v <f>, (4.77) 
rmpip — > — ah cos <p. (4.78) 

Our result shows that the transition point a c = l/(87r) corresponds to G c = —n/2. 
By making use of this equivalence, the Wilson loop ( |3.64| ) in the topological non- 
trivial sector of the four-dimensional U(l) gauge theory can be calculated in the 
two-dimensional massive Thirring model. The details will be given elsewhere. 

4.2 Lower and higher dimensional cases 

Using the equivalence between the TQFT part of U(1)d gauge theory and 0(2) D _ 2 
NLSM, we can study other dimensional cases. 

For D = 3, the equivalent 0(2) NLSM is one dimensional, 0(2)!. This is not a 
field theory model, but a quantum mechanical model of the plane rotor. There is no 
phase transition in this model. This implies that the three-dimensional U(l) gauge 
theory has only the confinement phase. This will be understood as the tunneling 
effect among classical vacua in the sense that the double- well anharmonic oscillator is 
related to the one-dimensional Ising model. It will be interesting to see agreement (or 
disagreement) on confinement mechanism between our approach and the Polyakov 
approach |2], [|. 

For D = 5, the 0(2) NLSM is three dimensional, 0(2) 3 . The three-dimensional 
0(2) NLSM has two phases on the lattice The phase transition is first order. 



This result is consistent with the mean field study of five-dimensional lattice U(l) 
gauge theory [Q. This theory has a finite non-zero critical coupling g c . In the strong 
coupling phase, quark confinement is expected to occur. However, the phase transition 
is first order. Therefore, on the lattice, it is impossible to take the continuum limit at 
this point. In view of this, the construction of continuum U(l) gauge theory from the 
lattice regularized theory is problematic in five dimensions and higher dimensions, 
unless the action is modified. 
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A Equivalence between Coulomb gas and sine-Gordon 
model 

The partition function can be rewritten as 

(ah) n 



Z sG (h) = I \d^}e- a J ****ffl>«»)) a 'g^- f d d xcos4>{x) 

J n=0 n - \-J 



d d x{e i4>[x) + e- 1 ^) 



2 n. 



ah) 2n 1 (In 



% (2n)\ 2*\nJ& 
(ah/2) 



n . 

J] / d d x i d d y i e i ^ ) - l,f,(n) 



2n n 



f [d0] e -«/^3(a^(»)) a ]T WW f[ f d d Xl d d yi e^-^\ (A.l 



Note that 

[[ e / ddx J( x )H^ 



[d<f>] exp 
1 



-f / d d x{d^{x)f + I d d .r.J [.,■)< 



exp 



— Jd d xJd d yJ(x)A(x iy )J(y) 



where A(x,y) is the massless scalar field propagator, 



A{x,y) := J 



(2ir) d p 2 



In particular, for J(x) = i Y,j QjS(x — 



3=1 



3 

exp 




x 



3)1 



Ei,* WkH x 3> x k)] for ?7 = 
for 3j 7^ 



(A.2) 



(A.3) 



(A.4) 



where the latter case is a result of invariance under constant translation of the field 
(p(z) — > </?(:r) + c. Using this result for qj = ±1, 



2n n 



( n! ) 2 i= i 



1=1"' ^ i=l 



{ah/2) 2n " 

- ^ A(xj -%)]}. 



n „ i 

E II / d d xd\exp { - —[J2l A ( x i ~ x i) + - y s )\ 

n=0 V 1 -) i=l J ZO£ i<j 



(A.5) 
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B Estimate of sine-Gordon partition function 



Note that the quantity 

g{xf - hg{x){-d 2 + hy l e {x) = e (x) _~ +h Q(x) (B.l) 
is positive, since (— <9 2 )/(— <9 2 + h) is a positive operator. Therefore, 

e -f(^/ s W-M I )(-a 2 + r 1 e W} ) g{x)= r M_ m{x) (B.2) 

is monotonically (rapidly) decreasing in G S}. Therefore, in the partition 

function, 

Z sG [r)} = £ ex P [ - IT ( 27r ) 2 / d d ^f - hg(x)(-d 2 + h)' 1 g(x)}],(B.3) 

{m(x)eZ;x€R d } 

the most dominant contribution comes from a set of configurations {|^(a;)|;x G S} 
which gives the smallest value for / d d x{g(x) 2 — hg(x)(—d 2 + h)^ 1 g(x)}. 

If the argument x of i](x) is outside the loop C{C = dS), i](x) = 0, while i](x) = 
2nq/g if x is inside the loop. For Z s g[0], {m(x) = 0} gives the most dominant 
contribution. Hence, we see 

Z sG [0] = 1 + ---. (B.4) 

For Z s g[i]}, the most dominant contribution comes from a set of integers {m(x)} 
whose value is the nearest to q/g where g(x) = 1-n\qj g — m(x)\. Since the integral 
part of q/g is absorbed in the shift of m(x), it is sufficient to consider the case 
< q/g < 1 without loss of generality. For a half-integer q/g, i.e., q/g = ±1/2, we 
see that {m(x) = 0, ±1} give the smallest value of 2n\q/g — m(x)\ for x G S. For 
< q/g < 1/2 (resp. 1/2 < q/g < 1), the most dominant contribution is given by 
{m(x) = 0} (resp. {m(x) = 1}). Thus, the rough estimate, for example, in the case 
of < q/g < 1/2 leads to 

Z sG [rj] = e-^Js^Mxf-hvW-dZ+h)-^)} + ..._ (B>5) 

A naive estimate of the ratio Z sG [r]]/Z sG [0] is given for r\ 2nZ by 

(WcM)sa=^^ = e-^Is^^-H-B^H x) + ^^ (B6) 
Z s gM 

Here (— (9 2 + h)~ 1 (x, y) is the massive scalar propagator with mass m ~ \fh and hence 
has an exponential damping factor e~ m \ x ~ y \. Therefore, the integral f s d d xrj(x)(— d 2 + 
h^rj^x) converges to a finite value even for large S. The term J s d 2 xr](x) 2 is pro- 
portional to the area of S. This term leads to the area decay of the Wilson loop 
expectation. 
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